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Abstract: We study perturbative and non-perturbative aspects of heavy-quark frag- 
mentation into hadrons, emphasizing the large-x region, where x is the energy fraction 
of the detected hadron. We first prove that when the moment index and the quark 
mass m get large simultaneously with the ratio NA/m fixed, the fragmentation function 
depends on this ratio alone. This opens up the way to formulate the non-perturbative 
contribution to the fragmentation function at large iV as a shape function of m(l — x) 
which is convoluted with the Sudakov-resummed perturbative result. We implement this 
resummation and the parametrization of the corresponding shape function using Dressed 
Gluon Exponentiation. The Sudakov exponent is calculated in a process independent 
way from a generalized splitting function which describes the emission probability of an 
off-shell gluon off a heavy quark. Non-perturbative corrections are parametrized based 
on the renormalon structure of the Sudakov exponent. They appear in moment space as 
an exponential factor, with a leading contribution scaling as NA/m and corrections of 
order {NA/mY and higher. Finally, we analyze in detail the case of 5-meson production 
in e+e~ collisions, confronting the theoretical predictions with LEP experimental data 
by fitting them in moment space. 



1 Introduction 



The fragmentation process, in which an energetic quark becomes a hadron, is one of the 
most interesting processes in the physics of the strong interaction, making an immediate 
hnk between the perturbative regime and confinement. The case of heavy (bottom 
and, possibly, charm) quark fragmentation is special in that the quark mass m, being 
significantly larger than the QCD scale A, allows one better theoretical control. In 
particular, the mass provides a physical infrared cutoff for collinear radiation making 
the fragmentation process accessible to perturbative methods. 

The methodology to describe processes involving fragmentation is based on factor- 
ization. In general "factorization" refers to the separation of the physical process into 
subprocesses which are mutually incoherent each depending on a separate scale. In 
practice this term is mainly used in the context of separating the perturbative and non- 
perturbative contributions p. The general procedure (see e.g. 0) is then to compute the 
cross section perturbatively and complement it by a convolution with a non-perturbative 
fragmentation function^ describing the softening of the heavy-quark momentum as it goes 
through the hadronization process. 

Experimentally, observables involving bottom fragmentation are particularly impor- 
tant. New experimental data have recently provoked much debate [SlIllE] about the 
accuracy of predictions for 5-meson hadroproduction obtained by convoluting the per- 
turbative cross section for bottom quarks with commonly used phenomenological mod- 
els |n|[7|. There exist many different viable implementations of a perturbative calculation 
for heavy-quark production using massive or massless quarks, different orders in a<(, soft- 
gluon resummation, Monte-Carlo simulations, etc. The phenomenology is particularly 
intricate because the parameters of the non-perturbative fragmentation function deter- 
mined by fits to data are very sensitive to the details of the perturbative description 
used. Consequently the application of these functions requires to use the same kind of 
perturbative description with which they were determined. 

On the theoretical side progress has been made at both the perturbative [Hlinilllll 
^lElEl and the non-perturbative P^ [T3] frontiers. Nevertheless, the most pressing 
questions have not been answered, namely 

• How should the non-perturbative fragmentation function be parametrized, in par- 
ticular in the region where the energy fraction x of the detected meson is large? 
This is especially important since the distribution peaks at large x. 

• Is this function really universal? e.g. can the parameters be fixed from data 
on bottom fragmentation in e~^e~ annihilation and used in various observables in 
hadron colliders? 

^Note a possible source of confusion in the terminology: the term "fragmentation function" is of- 
ten used both for the single-particle inclusive cross section in e^e~ annihilation and for the function 
describing the hadronization of a parton into an observed hadron. In this paper we use the latter. A 
fragmentation function is not an observable. When evaluating the cross section for a given physical 
process, the fragmentation function is convoluted with the proper coefficient function. 
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Answering these questions requires deep insight into the non-perturbative regime. 
However, a crucial step in addressing them is making a clear and sensible separation 
between the perturbative and the non-perturbative aspects and then dealing with both. 
Most previous work on the subject is sharply divided between a purely perturbative treat- 
ment, which concentrated on resumming the dominant perturbative corrections [HEIIlEl, 
and a purely non-perturbative approach, which disregarded the perturbative aspects. 
The latter was either phenomenological in essence [HllZl, or relied on general consider- 
ations such as the quark-mass expansion |l4j|, but avoided the large- a; region. None of 
these approaches has led to a satisfactory description of the fragmentation function at 
large x. This region is characterized by large perturbative and non-perturbative effects, 
so both aspects need to be addressed. 

Owing to soft gluon radiation the perturbative result at any given order diverges for 
X — > 1. This makes Sudakov resummation an essential ingredient in the description 
of heavy-quark fragmentation fl^.. While important, Sudakov resummation (see e.g. 
Fig. IHl below) does not lead, by itself, to a phenomenologically acceptable description 
of the cross section: it is much too hard and eventually becomes negative near x = 1. 
Non-perturbative corrections are required. 

It is well known that due to renormalons ^Hl the separation between the perturbative 
and non-perturbative regimes is arbitrary. This means |17 [ ni ^ lT? H I^ on the one hand 
that great care should be taken in controlling the accuracy of the resummed perturbative 
expansion whenever power-corrections are important but, on the other hand, that some 
insight on non-perturbative physics can be gained based on perturbation theory. Under 
the assumption that the renormalon contribution dominates, this information can be used 
to parametrize the non-perturbative corrections. Realizing this, Nason and Webber p!3] 
have used the structure of infrared renormalons to determine the parametric form of the 
leading power correction to the heavy-quark fragmentation function in the large- limit, 
where N is the moment (at large the A^-th moment is sensitive to x ~ 1 — 1/A^). 
They found that it should scale as NA/m. This is an important result, which we further 
establish and generalize here. 

Based on the quark-mass expansion of jHj and on a general representation of the 
fragmentation function at large A^ in terms of a non-local lightcone matrix element at a 
large lightcone separation (which follows directly from the definition P), we prove that 
in the limit where A^ and m become large simultaneously, the fragmentation function 
depends only on the ratio NA/m up to corrections of order 1/A^. This opens up the way 
to resum these non-perturbative corrections into a function of a single variable NA/m. 
The corresponding function in momentum fraction space depends on m(l — x). This 
function is convoluted with the resummed perturbative distribution, thereby introducing 
a shift and a deformation of the shape by non-perturbative effects. For this reason it 
is called a "shape function". The scaling property of the fragmentation function was 
known before (see j2] are refs. therein), but a formal proof based on the properties of 
the hadronic matrix element was never established. 

Stepping forward from a leading, additive power correction to a shape function, which 
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is a multiplicative correction in moment space, is essential for describing differential cross 
sections near a kinematic threshold pT|l22H23llTH] . Shape-function based phenomenology 
has been developed particularly in the context of event-shape distributions in e^e~ an- 
nihilation 121112211211^] where it has proven very useful. The same methodology can be 
applied to a large class of observables. 

Our approach to the problem of heavy-quark fragmentation is based on perturbation 
theory. Instead of parametrizing the differential distribution for the production of a 
heavy meson from a heavy quark, we start off with a perturbative calculation of this dis- 
tribution. The result is then modified by power corrections through a convolution with a 
shape function whose parametric form is deduced from the ambiguity in the perturbative 
result. A convolution of a perturbative (possibly Sudakov-resummed) heavy-quark dis- 
tribution with a phenomenological non-perturbative function has been commonly used 
in describing heavy-meson production. However, there are deep conceptual differences 
between our approach and the common practice: a) the resummation we perform, which 
focuses on the effect of the running coupling in the Sudakov exponent, is aimed at power 
accuracy, so that the separation between the perturbative and non-perturbative compo- 
nents is power-like; b) since the functional form of the non-perturbative component is 
constructed according to the ambiguity of the perturbative one, it should not be con- 
sidered as a phenomenological model; c) since we only rely on QCD, there is room for a 
systematic improvement of both the perturbative and the non-perturbative ingredients. 

We concentrate in this paper on the large-x region, neglecting power corrections which 
are suppressed by the energy of the quark, 0{a1) perturbative corrections that are not 
logarithmically enhanced, as well as power corrections in A/m which are not enhanced 
by the same power of N. Consequently, our description of the first few moments is of 
limited accuracy. 

Attempting to reach power accuracy at large x, our resummation program is based on 
Dressed Gluon Exponentiation (DGE). Calculating the Sudakov exponent as an integral 
over the running coupling (a renormalon sum) we refrain from making a truncation at 
some fixed logarithmic accuracy. The advantages of this methodology have already been 
discussed in detail in [THlliniEn] • However, some new features appear in the case of the 
fragmentation function owing to the low scale m/N. As we will see, the perturbative 
expansion of the Sudakov exponent reaches the minimal term already at the next-to- 
leading-log (NLL) and a rather large difference exists between truncation of the series 
at this order and a principal-value (PV) regularization of the Borel integral. The latter 
is advantageous as a basis for the parametrization of power corrections, being it free of 
Landau type singularities and having a smooth behavior at large A^. 

Although data analysis is performed here only in one case, namely the semi-inclusive 
cross section for i?-meson production in e~^e~ annihilation at LEP, we are guided by 
the hypothesis that the fragmentation function is universal, an idea that guided many 
phenomenological applications in recent years |2nil2Sll2Zll2HllS]- With universality in 
mind, we write our resummation formulae for the cross section in a fully factorized form, 
where subprocesses (such as the fragmentation function) depending on distinct external 
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Figure 1: The definition of the fragmentation function. The dashed hne represents a 
path-ordered exponential. 

scales appear as independent entities. This type of factorization is not restricted to the 
resummed perturbation theory, but instead it becomes a property of the full subprocess 
including the non-perturbative contributions. 

The paper is organized as follows. Section 2 recalls the standard definition of the 
fragmentation function and then deals with the heavy-quark fragmentation function in 
the large- limit on general grounds. Section 3 summarizes a process-independent, all- 
order calculation of the heavy-quark fragmentation function in the large-/9o limit. It is 
based on computing a generalized splitting function describing an off-shell gluon emission 
off a heavy quark. We then analyze the Borel singularity of the fragmentation function 
and construct an ansatz for the shape function. Section 4 specializes to the case of e^e~ 
annihilation. We begin by performing a renormalon calculation for the differential cross 
section of the entire process and then identifying the contributions which are enhanced 
near x = 1. Isolating the fragmentation subprocess we recover the result of Sec. 3 and 
interpret the remaining subprocesses which form the coefficient function. In Sec. 5 we 
discuss the phenomenological implications of our results and confront them with LEP 
data on bottom production in e^e~ annihilation. Section 6 summarizes our conclusions. 

2 Heavy quark fragmentation in the large-x limit 

The leading-twist quark fragmentation function, corresponding to the probability of a 
quark to fragment into a hadron with a longitudinal momentum fraction < a; < 1 (see 
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fig. Hj) , is defined^ by P , 

D{x;^i^) = -^[ — exp(zpi?//x)F(pi2/;/i2), (1) 

where F(piy] /i^) stands for the following non-local matrix element 

F{piy;fi') = E Tr{{O\0iy)\Hip^)+X){Hip,)+Xmom,^}, (2) 

renormalized at /i^, averaged over spin and color of the initial quark and summed over all 
final states containing a hadron H with momentum pi. We choose light-cone coordinates 
such that pi has a large "+" component, while y is a light-like vector in the "— " direction 
(?/+, y± = 0). Path-ordered exponential factors are required to make © gauge invariant. 
However, since we will eventually use the light-cone gauge y ■ A = where these factors 
become trivial, we will not write them explicitly. 

While this definition applies to fragmentation of both light and heavy quarks, there 
is a fundamental difference in our ability to use the theory in the two cases. In the case 
of light quarks, only the /i^ evolution of the matrix element is calculable, whereas for 
heavy quarks with mass m ^ A the matrix element itself and thus m^;/x^) can 
be calculated perturbatively, while non-perturbative effects can be treated systemati- 
cally as corrections. The matrix element can be calculated since provides a natural 
cutoff for collinear gluon radiation. It also provides a natural scale for renormalization 
of the operator /i^ = m^, serving as a factorization scale in a generic hard process. 
In the following we will use simplified notation D{x,m'^; fi^ = m^) = D{x,m'^) and 
F{piy, m^; /i^ = m^) = F{piy, rn?). 

By matching the general definition given above onto the heavy-quark effective theory, 
Jaffe and Randall ^1] have shown that in the large m limit the non-local matrix element 
in Q can be written in terms of a function of the product {piy) A/m up to corrections 
which are suppressed by A/m: 

^^^py"^^^ i-Wiy) [j^iPiy A/m) + 0(A/m)] . (3) 

Physically A appears as the energy of the light quark (including the binding energy) in 
a meson, corresponding to M — m, where M and m are the heavy meson and quark 
masses, respectively. This particular dependence of F{piy,rn?), which follows directly 
from QCD, has a phenomenological significance: it constrains possible models for the 
parametrization of the heavy-quark fragmentation. Due to the difficulty to distinguish 



^It should be noted that, in contrast with F(j)iy\ jjt^) is Lorentz-frame dependent. Eqs. Q 

and (O refer to a Lorentz frame in which the decaying quark has a vanishing transverse momentum, 
in accordance with the Sudakov parametrization (|17|l we will use below. It is for this reason that 
Eq. Q differs by a factor a;'^^^ (where d is the space-time dimension) from the standard definition ^ 
corresponding to a frame where the hadron (pi) has a vanishing transverse momentum. 
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between perturbative and non-perturbative corrections at large x, Jaffe and Randall 
exclusively concentrate in their analysis on the first few Mellin moments of the D{x, m?), 

b(N,m^) = dxx^-^D{x,m^). (4) 
Jo 

As explained in the introduction, the behavior of the fragmentation function in the 
large-x limit is particularly interesting both from a purely theoretical point of view and 
for practical applications. As we show in the following sections, if properly resummed, 
perturbation theory does provide a basis for the analysis of heavy-quark fragmentation 
at large x, so long as m(l — x) ^ A. Of course, power corrections can be ignored only 
if m(l — x) ^ A, so in practice the perturbative calculation must be supplemented by 
non-perturbative corrections. 

While phenomenological parameters are definitely needed for a precise description 
of (^, it is important to constrain theoretically its functional form as much as possible. In 
the rest of this section we show how the large-x behaviour of -D(x, /x^) can be constrained 
based on the Jaffe and Randall result (jHl). 

First, let us rewrite the definition of the moments Q for ^ 1 as 

D{N,m^)c^ dxexp{-N{l~x))D [x,m ), [0) 

Jl-No/N 

where we used the fact that the integral is dominated by the x — > 1 region to expand 
In X ~ — (1— x) and to modify the lower integration bound to 1—Nq/N with 1 <^ A"o <^ A^. 
Corrections of order 1/A^ are neglected here. A possible choice of A"o is This 
guarantees that the region x ~ 1 — 1/A^ dominating the A^-th moment is fully 
included in 1^. 

Substituting into © and changing the order of integration we obtain, 
D{N, m^) ~ — / d{v\y) \ — exp (-A^(l - x) + %V\ylx) 



271 J-oo Piy Jl-No/N X 

77-/ d{piy) ^^P^y^^ ) exp{ipiy) [ d^ exp {{- N + ipiy) ^) , (Q) 

111 J-oo P-IV Jo 



where in the second line we used the large x approximation and changed the integration 
variable to ^ = 1 — x. Performing the ^ integral we get 

[exp (-ATo + ^p,yNo/N) - 1] 'J^^^. 

ipiy - N 

In the first factor we can safely take the limit A^ — > 00 and then A'o — >■ 00 removing 
the Aq dependence completely. The result is: 

D{N, m") r d{p,y) EiMlJlfl L,^ exp {tp,y) , (7) 
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This last integral can be performed by closing a contour in the complex piy plane, 
assuming that the matrix element F{piy, im?) /piy itself has no singularities. This integral 
should be done with some care: the phase of F{piy,m'^)/piy dominates that of the 
exponential factor, so the contour should be closed in the lower half plane. 

To verify that the phase of the matrix element dominates, consider the inverse Fourier 
transform to (jT]): 



F{piy,m? 
Piy 



1 dx 



X 



D(x,m^) exp{—ipiy / x) 



(8) 



where we used the fact that D(x,m^) has support only between and 1. Now, since 
D{x^m?) is a positive definite function, Eq. (jSJ represents a weighted average of the 
phase factor exp(— zpiy/a;). This phase is opposite in sign and larger than ipiy for 
any < x < 1. 

Having establish this, the integral (jZj) can be evaluated by closing the integration 
contour in the lower half plane (the integral along the contour at infinity vanishes), 
picking up the pole at piy = —iN. The result is the following remarkable relation 



D{N, 



F{piy, 



exp {ipiy) 



+ 0{l/N). 



(9) 



piy=-iN 



Eq. Q associates the behavior of the fragmentation function at large N and that of 
the non-local matrix element (j21), analytically continued in the complex piy plane and 
evaluated at asymptotically large light-cone separation piy = —iN. We recall that 
similar asymptotic relations exist in deep inelastic scattering between light-cone matrix 
elements and moments of twist two jHS] and twist four [SB] distributions. 
Using now the result by Jaffe and Randall Q we deduce that 



piy= 



-iN 



+ 0{1/N), 



(10) 



namely, that in the limit where m and N get large simultaneously D{N,m?) becomes, 
asymptotically, a function of a single parameter NA/m, up to corrections suppressed by 
or by A/m. It should be noted that ^ holds at large m for any piy (not necessarily 
large), and likewise, holds at large for general m. However, from the two equations 
together it follows that there is one, specific way to take the simultaneous limit getting 
a function of a single argument, namely with the ratio N/m fixed. 

We recall that Ref. ,14] used Eq. to derive a scaling law for the fragmentation 
function in x space. They also performed moment space analysis but concluded that 
their result is only applicable to the first few moments, due to the divergence of the 
expansion in A/m. This problem is avoided, however, upon taking the limit where m 
and get large simultaneously. 

The appearance of the scale m/N, or in x space, m(l — x) is not surprising. As 
we discuss below, this scale emerges naturally in perturbation theory: m(l — x) is the 
transverse momentum at which soft-gluon radiation from the heavy quark peaks. It is 
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directly related to fact that the "dead cone" ^ angle is proportional to the ratio between 
the quark mass and its energy (see Sec. 3.2.1). 

Let us also recall that the scaling of the fragmentation function at large x as m{l — x) 
was discussed in the literature in the past. A particularly intuitive argument in favour 
of this scaling law was provided in It was shown there that this behavior follows 
from simple kinematic considerations if one assumes that the differential cross section 
depends on the mass of the heavy quark and on its energy only through the ratio, i.e. 
through its velocity. 

The presence of a correction 0{NA/m) was also predicted by ^Hj based on a renor- 
malon calculation in the large- limit. In what follows we will adopt the assumption 
of ^3] concerning renormalon dominance at large A^. We will be using the renormalon 
structure of the Sudakov exponent P^ll^ to construct an ansatz for T{piyK/m). As we 
explain in Sec. 3.2, by excluding an additional logarithmic dependence on m, Eq. (fTOl) 
has important consequences concerning the nature of the renormalon singularities in the 
Sudakov exponent beyond the large-/5o hmit. 

3 Heavy- quark fragmentation function — a process- 
independent renormalon calculation 

3.1 An ofF-shell splitting function 

Consider gluon emission off an outgoing heavy quark in a generic process. Assume that 
the quark is on-shell after the emission p\ = m^, while the gluon has some time-like 
virtuality k"^. 



Figure 2: Single gluon emission off a heavy quark involved in a generic hard process M-q. 
The amplitude for the emission of a single gluon off this quark (Fig. ^ is 





(11) 
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where is a colour matrix, e^* is the gluon polarization vector and Aio represents the 
rest of the process. 

We wish to compute radiative corrections which are associated with the singularity 
of the quark propagator + k)"^ — m^] = l/{2pik + k"^). These corrections become 

dominant for a class of observables, including the fragmentation function. Choosing the 
lightcone axial gauge n ■ A = 0, 

d,. - E ^t*4 = -9.^ + (12) 

A 

where n is some^ lightlike vector (n^ = 0), interference terms are suppressed and the 
entire singularity is contained in the amphtude (fTTj) squared, 

E -^-^^ = ,0 ?^!^'/2^2 ^M-Tr {(i/i + ^ + m)YiA + m)Yi]^i + ^ + m)A^oA^o} • (13) 

Here a sum over the quark and gluon polarizations was taken. 
Proceeding with the evaluation of (fT^ we write 

where ri 2 correspond to the two terms in the propagator ()12j). 
n = -Tr + ^ + m)7^(2/i + m)-ff'{^^ + ^ + m)MoMo} 

rs = Tr {(^1 + ^ + m) [?/(]/i + m)^ + ^(J^ + m)?/] + ^ + m)MQMo} . (15) 
After some algebra we get 

n = -4m [m^ + pik + k^] Tt{MoMo) - 2 [2m'^ + k^] Tt{]^^MoMo) 
+4 [-m^ + Pik] Tri^MoMo) 
2{2'D^ k + k'^^ 

ra = ^^^-^7 {m[2pin + kn] Tr(A^oA^o) + ['^pm + kn] Tri^^MoMo) 

+ PinTr(^7Wo-Mo) -piA;Tr(7/7WoA^o)}- (16) 
Next we introduce the following Sudakov parametrization (adopting the notation of [T^). 



■^The vector n is usually set parallel to the direction of the other incoming or outgoing parton in Aloj 
however this is not essential. 
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where = rn? and k±_ is orthogonal to both n and p, with k^^k^ 
variables, 



-k"^, . In these 



n = -Am [m^ +pik + Tr(A^o^o) - 2 pm^ + A;2 _ (i _ 2)(2piA; + P)] Tr(2/A^oA<o) 
-2 (2piA; + Tr(^^A<oA<o) + ^^^'^^ [-2(2 - z)m' + z(2pifc + A;^)] Tr(7/A<oA<o) 

2z 



j9n 



r2 = 2{2pik + k^) \m 



l + z 
I- z 



Tr(A^o-Mo) 



1 - z 



1 — Tn — 

Tiim^Mo) + — Ti{4MoMo) 

1 — z pn 



(18) 



and finally, 



(2piA; + k^y 



2m 



2z 
1-z 



{2pik + k'^) - (2m2 + fc^ 



Tr(;Wo;Wc 



+ 2 



1 + ^2 



+ 



1 - z 

2pik + k^ 
pn 



{2pik + A;2) - (2m2 + A;^ 



Tr(j/^o-Mc 



2z 
1-z 



(2piA; + k^) Tr(I^^A^o-Mo) 



[z(2piA; + A;2) - 2(1 - z)m^] Tr(?/A^o^o) 



(19) 



The quark propagator has the following singularity: 

1 1 z{l-z) 



[pi + kY — m^ 2pik + /c^ zk"^ + k\ + [1 — z^m^ 



(20) 



Therefore, the relevant limit^ is the one in which m^, k\, k"^, all become small simul- 



taneously, whereas the ratios between them (which depend on the quark longitudinal 
momentum fraction z = (pin) / {pin + kn)) are fixed. 

In this limit the dominant contribution to the squared matrix element ()19|1 is the one 
proportional to Tr(]/7Wo-^o)- All the other terms are suppressed by one of the small 
parameters. We thus proceed with the following approximation. 



{2pik + A;2)2 



l + z' 
l-z 



{2pik + k'^) - {2m^ + k'^ 



Tr(2/>[o-Mo), (21) 



which is factorized into a product of the squared matrix element of the non-radiative 
process M.qM.q) and a splitting function describing the probability of a single gluon 
emission. This splitting function corresponds to the emission of an off-shell gluon off a 
heavy quark, and thus it generalizes both [THj where the gluon is on-shell and jSH] where 
the quark is massless. 



*This is a generalization of the quasi-collinear limit ^21 to the case of an off-shell gluon. 
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3.2 The quark fragmentation function in the large-/5o Umit 

Based on the general definition and the off-shell splitting function derived above, 
we can now calculate the quark fragmentation function in a process independent way. 
The calculation is performed with a single dressed gluon, so it is exact at order ag and 
contains the leading contribution in the large-/3o limit at higher orders. Using DGE, this 
result will be used to compute the Sudakov exponent. 



3.2.1 NLO calculation 



Let us begin by extracting the Born-level result. Saturating the state \H{pi)+X) 

by the outgoing on-shell quark with momentum pi alone, and using the definition (jSJ 
one finds that 

FiPiy) = ipiv) exp{-ipiy). (22) 
The ?/_ integration in can be readily performed. 



11/"°° ■ f 

-- / d{p,y)e'M 

X 2lT J-oo 



5(1 - a;) 



identifying the longitudinal momentum fraction X as 1 . 

The NLO {0{as)) calculation proceeds in a similar way. The state \H{pi) + X) 
is saturated by the outgoing on-shell quark pi and off-shell gluon k. Calculating the 
squared matrix element (j21) we get (PT|l with TVIqA^o = exp{—ipiy/z), so 



F{Piy) 



Piy ^j,,yi^ Stt Cpas 
z {2pik + k^f 



l + z^ 
1 - z 



{2pik + - (2m2 + A;^ 



(23) 



where d(j) is the gluon phase space with the longitudinal momentum of the "detected" 
quark pin fixed. We see that the state \H{pi) + X) is effectively replaced by a single 
outgoing on-shell quark with momentum pi/ z, times some factors which are incorporated 
into the splitting function. 

Let us perform this integration with a fixed gluon virtuality. The phase-space integral 

is 



d^k 
(27r)3 



5{2k+k^ -kl-k^ 



167r2 



dz 



z(l 



dk' 



(24) 



The contribution of a single real gluon emission to the quark fragmentation function is 
therefore 



D'''^\x,k^/Tn^) 
Cpas 



2tt 



dki 



2n x{l — x) 
1+x^ 1 



dki 



1 + x^ 1 {2m^ + P) 

1 - X 2pik + k^ ~ {2pik + A;2)2 

(2m2 + e) 



(25) 



1 — X xk"^ -|- /c^ + (1 — x)'^7ii? 



— xil — x) 



(xk^ + k'^, +(l-x)2m2)2 



where we substituted the expression for the propagator (pUj). 
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Next, performing the integration over k'j_ from k'j_ = (ignoring terms from the upper 
integration hmit where the propagator is non-singular), we get 



CfQs 
27r 



— — In (x ik"^ Im?) + (1 — X 
1 — X ^ 

, 2 + (fcVm^) 



xil 



X [k'^ /w?) + (1 — xY 



(26) 



Setting fc^ — i> in this expression we recover Eq. (58) of [TH], as expected. 

Let us return now to the non-integrated form of the real-emission matrix element ()25j] 
and examine the distribution of the radiation in some Lorentz frame where the longitu- 
dinal momentum of the quark pn = pj^ = y^/v^ is large: 2> fc^, k^^rn? {\f(p can 
be, for example, the centre-of-mass energy in e^e~ annihilation). The energies of the 
outgoing quark and gluon are 



q X + k + 



pi 



Eu 



2xvV 2 
g2(l - xf + k'^ + k] 



.2 



2{l-x)^ 



■x + o[l/^^ 



(27) 



so in this frame the longitudinal momentum fraction coincides with the energy fraction 
and the emission angle is 



sin^^ 



q^x^{l — X) 



+ 0{l/q' 



(28) 



The leading (log-enhanced) contribution at 0{a, 
1 -x < 1, 

dD"'''' as sin^6' 



is obtained by taking A;^ = and 



dx d sin^ 6 



TT 



1 



(29) 



x)(sin^6' + 4m2/g2)2- 

This result implies that the radiation peaks at angles 9 ~ 2m/ and that it is depleted 
in the forward direction (^^ = 0). The authors of Ref. [H] called this depletion the "dead 
cone", contrasting it with the enhanced radiation for 9 — > in the case of a massless 
quark. In a boost-invariant formulation the radiation peak appears at \k^\ ~ m(l — x). 
We have already seen in Sec. 2 that m(l — x) ^ m/N is the relevant scale at large x. It 
will emerge again as the typical gluon virtuality (/c^) and eventually, as the natural scale 
for the coupling in the framework of the renormalon calculation that follows. 



3.2.2 All-order calculation in the large-/5o limit 

The 0{as) result for an off-shell gluon emission derived in the previous section can be 
promoted to a single dressed gluon (SDG) renormalon sum by integrating over the gluon 
virtuality k"^, with the coupling being renormalized at k"^ (see fig. E]). It is convenient 
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y 




Figure 3: The diagram contributing to D{N, im?) in the hght-cone axial gauge A-y = Q 
and in the large Nf limit. 



to define a factorization scheme and scale invariant quantity, by taking the logarithmic 
derivative with respect to m^. Contrary to the fragmentation function itself, its logarith- 
mic derivative is free of ultraviolet divergence. The result, written in moment space Q 
scheme invariant Borel transform jSD], is 



du T{u) (mVA^) " f dx (x^~^ - l) Bd{u,x), (30) 



d In m? Po Jo 

where x^~^ has been replaced by x^~^ — 1 to account for virtual corrections, and 

„ , , -MSinVTM ^„ 2« /o-.\ 

Bd{u, x) = e''" m^" (31) 

2n 

1 + ^' , f n2, 2^ , ^2^ , ^ 2 + (P/m^) 



X [k"^ /m?) + (1 — x)^ 



Note that in m? can be scaled out, so B£){u,x) is independent of any scale and the 
dependence of (ilnDsoc(A^, m^)/(ilnm^ on appears only through {m? / h?) Here 
the factor (— m) is associated with the logarithmic derivative, the sine appears upon 
taking the time-like discontinuity of the dressed gluon propagator, and T{u) is defined 
as the Laplace transform of the coupling 

/■oo , , _^ 

A{k^)= (3^as{k^)/'K = duT{u){k^/X^) . (32) 

T{u) depends only on the coefficients of the renormalization group equation of Os/tt. For 
a one-loop running coupling T{u) = 1 and in the two-loop case 

T{u) = (u5)"^ exp(-n(5)/r(l + u6), (33) 
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where 6 = Pi/Pq, with Pq = vrfoo and Pi = ti%i as defined in Eq. (28) of [TSI. Defining A 
in the MS scheme, 

c = 5/3 + 0(l//5o). (34) 
Performing the integration over the gluon virtuahty in (|HT|) we get 



Bd{u, a;) = — e'^ 



X 



Here we used integrals of the form 

dy 1 



x) (l + u) 



(35) 



TT 



y 



1 



x] 



sm Tcu 



x) 



-2{l+u) 



with y = xk'^/m?. Note that this sinvr-u factor in the denominator cancels against the 
one in the numerator of (j3T|l . so there are no renormalon singularities in (j35|) . There is 
however, a convergence constraint for the Borel integral (jHn|) at u — > oo. In the case of 
one-loop running coupling {T{u) = 1) the constraint is: 

m^(l — X f/x>A^exp{c). (36) 

As we will see below, this singular behavior for x — > 1 translates into renormalons upon 
performing the x integration. 

Note also the appearance of the factor (1 — n) in front of the x — > 1 singular term 
in ()35p. This refiects a relation between the two terms in the squared matrix element ()2ip. 
Let us examine the two corresponding terms in the square brackets in ()3ip : the first is 
the source of the 1 and the second, which can be represented as a logarithmic derivative 
of the first with respect to k"^, yields the —u. As we will see, in moment space this 
structure leads to the absence of a renormalon singularity at m = 1, while renormalons 
do appear at all other integer and half integer values. 

Let us perform now the integration over x in Eq. (j^O)). We obtain 



Po Jo 



du T{u) (mVA^) "5^(M,iV), 



(37) 



with 



dx {x 



N-l 



-e" 



[l-u)T{-2u) 



- Ij Bn{u,x) = 

' T{N + l + u) _ r(2 + M) 
r(A^ + 1 



u 



2u) 



T{2-u) J 

' r(A^ + n) 
r(iV + 2-n) 



r(i + u) 



r(3 



(38) 



As anticipated, Bi){u, N) has infrared renormalon singularities (simple poles) at all pos- 
itive integer and half integer values of u, with the single exception of m = 1. 
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3.2.3 Ultraviolet subtraction for the fragmentation function 

In order to obtain the fragmentation function itself Eq. p7jl needs to be integrated. The 
integration involves a subtraction of the ultraviolet divergent contribution, as standardly 
done in collinear factorization. The result takes the form: 



Cf 

Jo 



: 1 + 

duT{u) (mVA^ 



u 



-u B^{u,N) 



u 



(39) 



where is a factorization scale and B^{u, N) is the Altarelli-Parisi evolution kernel (the 
subscript E stands for Evolution) , 



7n(iV) 



n=0 



n! 



(40) 

by canceling the 



This subtraction guarantees the existence of the integral near u 
1/u singularity. 

We stress that this subtraction is required only because we calculate the fragmenta- 
tion process alone. In any infrared and collinear safe observable the coefficient function 
will compensate the 1/u singularity. This is demonstrated in Sec. 4 in the case of the 
differential cross section in e~^e~ annihilation. Through the evolution term the frag- 
mentation function becomes factorization-scheme and scale dependent although the full 
result for the cross section is not. 

The leading order coefficient in ()40|) is renormalization-scheme invariant, and it equals 
to the u = limit of ()38|) . ensuring the cancellation of the 1/u singularity in 



7o(iV) = S,{N) 



3 1 

4 + 2 



1 

N 



where ^^(A^) = J2j=i 1/j , so that ^i(A^) = ^(A^ + 1) + 7e- In the MS factorization 
scheme 7n(-^) are known to two-loop order (n = 1) in full 
large-/9o limit jHS]. The latter is given by^ 



and to all orders in the 



+ 



n=0 



'A^ + 2A-1 



^(A^ + A) -^(1 + A) 

1 {1 + Ay 



(41) 



1 



with 



A 



1 + A N + A 2 + A N + l + A^ 
sirnrA T{A + 2A) 



+ O{l/Po) 



TT 6T{2 + Ay 



+ O{l/(3o). 



(42) 



^Note that using the scheme invariant Borel transform beyond the large- /So hmit, the relation between 
the coefficients of Bj^{u,N) and those of E{N,A) becomes non trivial (see e.g. |2()|^. 
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A is the all-order large-/3o expression in MS for the gluon bremsstrahlung effective 
charge [32] or the cusp anomalous dimension PHIH^j lHU]. and A = PQas/ir is the large-/9o 
coupling in MS. 

Using eqs. (jH^ through pT|l it is straightforward to extract the large-/3o coefficients 
of the heavy-quark fragmentation function D{N, m?) in the MS factorization scheme to 
all orders, for example 



Sl{N)+ 

1 1 1 1 „ 

A{m^) + 



^ ^ +l')5i(iV) + l-52(iV)+ ^ ^ 



AT A^+i ' y ^ ' ' 2(iV + l) 2N (A^ + l)2 

' --Sl{N)+(- i -\ Sl{N) 

3 ^ VA^ iV+1 3/ ^ 




5^,,.. 1^,,.. /Ill 1 1 \ 2 173 19 1 19 1 



6 ^ ^ 6 ^ ^ V4 6A^ + 1 QNJ 96 36 36 + 1 



+ — + 



5 19 1 11 



12Ar2 4(A^+1)2 4A^3 4 (AT +1)3 



A^m'f + 0(A(m2)') . (43) 



The order as result is known. It was extracted by Mele and Nason [Hj, starting from 
the single inclusive cross section for heavy-quark production in e^e~ annihilation and 
subtracting the coefficient function for the massless case. More recently, this term has 
been computed in a process-independent way in [T2]. Our calculation extends it to all 
orders in the large- /5o limit. 

In addition to the resummation of large perturbative corrections, the all-order large- 
j3o result can be used to extract some non-perturbative information on the fragmentation 
process. As we saw, in the large-/3o limit renormalons appear always as simple poles, 
and are exclusively associated with the integration over x. It should be noted that the 
subtracted term (Eq. (jlHl)) has no infrared renormalons in factorization schemes such as 
MS. 

In ()38|1 the first singularity is located at m = 1/2 and it corresponds to 0{A/m) 
corrections. The second singularity is at m = 3/2 {0{A^ /m^) corrections) and higher 
singularities appear at all integers and half integers on the positive Borel axis. It is 
particularly interesting to note the absence of the renormalon at m = 1. As explained 
above this structure can be traced back to the relation between the single and the double 
pole terms in 2pik + k'^ in the squared matrix element ()21|) . 

Both the presence on the renormalon singularity at m = 1/2 and the absence of the 
one at M = 1 are consistent with the findings of Nason and Webber which consid- 
ered the large-A^ limit of the fragmentation function in e^e~ annihilation. We confirm 
these conclusions in a more general, process independent context (and for a generic A^), 
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and derive an expression for the Borel function allowing to combine resummation with 
parametrization of power corrections. It is not known, and certainly deserves further 
investigation, whether a renormalon at u = 1 does appear in the full theory. 



3.2.4 The large- a: region and the Sudakov exponent in perturbation theory 
and beyond 

At large x multiple emission plays an important role. The SDG result can be readily 
used to derive the exponentiation kernel in the approximation of independent emission. 
Extracting the log-enhanced terms from ()38|) and ()40|) one gets 



DGE Jo 

where ii was set equal to m. 



du 



+ InN 



u u 



T{u) (mVA^) " (44) 



A^) = - e^" (1 - u) V{-2u) [A^'" " ij + 0{1/(3q), (45) 

and B_a{u) = 1 + + . . . is the Borel transform of the cusp anomalous dimension ()42|). 

Note that we write here a perturbative expansion for \n.D{N, m^). This is not just a 
formal manipulation which holds to order Og, but rather a statement about the structure 
of higher-order terms. The logarithmically enhanced terms, as opposed to the expressions 
in (jnZj), and exponentiate: they can be written to all orders as 

{oo n+1 ^ 
n=0 k=l } 

Although (jlH) was derived strictly in the large Pq limit, i.e. leading-log (LL) accuracy, 
it can be generalized in a straightforward way to comply with NLL accuracy in the full 
theory [S^18ll29j by using the two-loop running couphng and replacing c of Eq. 
in the exponent e^" and in -B^(u) by 

c ^ 02 = 5/3 + (1/3 - 7rVl2)CA//?o. (46) 

It should be noted, though, that the functional form of the Borel transform is not known 
beyond the large-/3o limit, so the replacement of c by 02 in ()45j) is not proven. As 
discussed in a similar context in [201 the same logarithmic accuracy can be achieved in 
different ways, e.g. by replacing the factor e^^/^-*" in ()45|) by e^^/'^''"(l + u{a2 — 5/3)). We 
will assume here the simple replacement ()46|) . 

The renormalon structure of the Sudakov exponent, is similar to that of the full 
SDG result (jHH|l : the leading renormalon is at m = 1/2 and subleading ones appear at 
all integer and half integers m > 3/2. In accordance with the general property of the 
fragmentation function (Sec. 2) these corrections show up in on the scale m/N. The 
appearance of these powers corrections in the exponent leads to factorization, in moment 
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space, of the non-perturbative contribution to the fragmentation function. Moreover, it 
also imphes that these corrections exponentiate together with the perturbative logs. 

Based on the result of Sec. 2 we can further deduce that in the full theory, as in (jSJ, 
renormalons in the Sudakov exponent appear as single poles, involving no anomalous 
dimension. In case on non- vanishing anomalous dimensions 7^, power corrections are 
modified by logarithms, (/iiV/m)*(ln/i/m)'''*/'^^^'') where /i is a cutoff. This behavior 
contradicts Eq. ()10|) and therefore the anomalous dimensions must vanish. 

We conclude that non-perturbative corrections to the heavy-quark fragmentation 
process appear as 

D{N, m^) = Dp^{N, m^) D^^{NA/m), (47) 

where Dpx(A^, m^) is obtained by exponentiating the (regularized) renormalons sum ()44p 
and 

D^P iVA/m = exp -cui-f u^-f ( " +■■■ , 48 

which is characterized by the absence of the second power of (NA/m). Here u;„ are 
left as free parameters. These parameters are defined only within a given regularization 
prescription for the renormalons, corresponding to a definite separation between the 
perturbative and non-perturbative components in D{N,m?'). A physically meaningful 
separation would be the one based on a momentum cutoff, or, alternatively, a principle- 
value regularization of the Borel sum, which can be readily related to a cutoff |17j . 
Having performed the renormalon sum in DpT^[N ,m?), the parameters in D^p{N A/m) 
can be fixed by comparison with experimental data. 

Going over to x space, the product in (jTTj) becomes a convolution with a shape 
function depending on m(l — x). As discussed in the context of event-shape distributions 
in e+e~ annihilation p71l211l2Hll8j . this convolution generates a shift of the perturbative 
distribution proportional to uJi and a deformation of its shape by c^a and on. 

It should be stressed that both the resummation formula and the non-perturbative 
corrections in (pSj) are expected to dominate only for ^ 1. The first few moments are 
surely affected by perturbative and non-perturbative contributions on the scale m (as 
opposed to m/N), which do not necessarily exponentiate and were neglected here. The 
phenomenological implications of and (jlH|) will be discussed in more detail in the 
context of e~^e~ annihilation. 

4 Heavy meson production in e+e~ annihilation 
4.1 Calculation of the Sudakov exponent 

Let us consider now the specific process of heavy-flavor production in e~^e~ annihilation. 
To be concrete, let us concentrate on the vector current contribution to the process 
e~*"e^ — > QQg where the momentum of the virtual photon is q (the centre-of-mass 
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energy is and the quarks have momenta pi and p2, energy fractions x = 2piq/q^ 

and X = 2p2q/q'^ and masses pf = pi = rn? = pq'^/4:. 

As in ^H] we begin with the exact matrix element with an off-shell gluon, /c^ 
The result for the cross section is, 



1 da 
a dx dx 



27r 7F 
_e + p/2_ 



P 



(x + e + p/2)2 + {x + e + p/2f - 2p(l + e + p/2) 



(l + p/2)(l-x)(l-x) 



p/2 



x)2 



(49) 



It is straightforward to verify that when taking the quasi-collinear limit this expres- 
sion reduces to the approximate one we calculated in the previous section in a process- 
independent way. Note first that the propagator which becomes singular in this limit is 
inversely proportional to 2pik + k"^ = q"^ — 2p2q = q'^{l — x), where we used momentum 



= q — Pi — P2- As follows from Eq. (j^UI) . the quasi-collinear limit is the 
X, e and p (and k\/q^) all become small but the ratios between them 
are fixed. In this limit the leading contribution to P9|l takes the form: 



conservation: k 
one in which 1 - 



1 da 
a dx dx 



2ti 



l+x^ 



p/2 + e 



1 — X [1 — x) 



[i-xy 



(50) 



which fully agrees^ with 1)211) . 

Owing to the non-vanishing quark mass and to its inclusive nature, the differen- 
tial cross section da/dx is an infrared and collinear safe quantity. It can therefore be 
calculated in perturbation theory to any order, with no need to perform factorization. 

Here we perform such a calculation having in mind the experimentally interesting 
scenario with p ^ 1 as well as 1 — x ^ 1 (or N ^ 1). These conditions are easily 
met for charm and bottom production at LEPl energies, i.e. g ~ 91 GeV. In this limit 
logarithmically enhanced terms of either p or N need to be resummed. At the same 
time, in order to address the issue of power corrections, we must perform renormalon 
resummation. This calls for DGE. 



^To verify the consistency of the two expressions one should write p2 in terms of the Sudakov 
parameters we used in the previous section: 



P2 = bp" + a{2pn)n" , 



1-62 



b {2pnY ' 

where a and b were fixed by the conditions p^ 



- p/2 -^x^-p 
2{l-x) + p/2 ' 

^ and (fc + pi + p2)^ 



(51) 

Next, evaluat- 



ing Tiij^MoMo) ~ 8pp2 = 4to2(1 + b'^)/b, one finds that in the limit considered a ~ q^/{2pnf-, 
Tr(2/Alo7\/(o) = 4(7^ and x ~ z, where terms linear in one of the small parameters were neglected. Fi- 
nally, making these substitutions in 1)21(1 . writing the phase-space measure as d(t> = / (IBtt^) dz dx and 
dividing by the flux factor 1/(4(7^)^ one find that the approximate result of the previous section indeed 
coincides with ((5?1)l . 
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In this section we perform the calculation by integrating the full off-shell gluon matrix 
element fl49|) over the exact phase space. This calculation follows the steps of the first 
derivation of the DGE result for the thrust distribution pH], based on the characteristic 
function of 38/. It proceeds by identifying log-enhanced terms and deriving a Borel 
representation which generates these terms to all orders in the large-/3o limit. Finally, 
this Borel sum is used as the exponentiation kernel. 

First we integrate ()49|) over x, getting a characteristic function. The limits of phase 
space (see jTH]) are: 

(2-x)(l-x-e + p/2)±S 



X 



max / mill 



2(1 



x] 



p/2 



x^ 



pj{l-x 



ep. 



(52) 



As Fig. m shows, for a given p the phase-space shrinks as e increases. Note that even 
with a vanishing gluon virtuality the phase-space boundary does not reach the line 
X = 1 (with the exception of a point sX x = 1), where the matrix element is singular. 
Consequently the differential cross section da/dx is infrared and coUinear safe. The 




Figure 4: QQg phase-space boundaries in the x-x plane, with varying gluon virtualities 
e = (1 - ^)2i/io, for i = through 9, with p = 0.04. 



resulting characteristic function is: 

1 da 

-— e,x,p = 
a dx 



(53) 



with 



-(p + 2e) 

4 - 4^ + 2^2 ^ 8e - 4ee - 2p^ + 4e2 - 



24e 



2^(2 + p) 

p + 2e 2(p2 + 48e) 



In 



+ 



'2e + 2e-2e' + 2ee-pe + 2S 
,2e + 2e-2e2 + 2ee-pe-2H 

2(-4 + p)(-p + 4e)\ „ 



,e(2 + p)p2 pe (2 + p)p2(4e + p) (2 + p)p(4e + p)2 



(54) 



where ^ = 1 — x. These three terms are distinguished by the singularity of the corre- 
sponding term in the matrix element at x — > 1: T\ originates in a double pole, T2 in a 
single pole and JF3 in non-singular terms. 

Next we perform the integral over the gluon virtuality e in the range < e < emax, 
where the upper integration limit, 



1 

4 



v/;^- Jp + 4(i-x) 



(55) 



is the point where S = and x 
sin TTM 



J. The Borel function is obtained by computing: 



!?(«, X, p) 



e™/(M,a:,p), I{u,x,p)=u / liee "J'(e, x,p). (56) 



The exact integration is complicated. The full result can only be written in terms of 
hypergeometric functions. However, in the following we will not need the full result, as 
we are interested in the logarithmically enhanced terms alone. Assuming the hierarchy 
p <C; ^ -C 1 the leading terms in the Borel function are rather simple. The integrals 
corresponding to the three terms in the characteristic functions are the following: 



smvru 
1 



(p/4)-", 



u 



TT 



SmTTM 
1 



r'-'"(p/4)-" 



1 



u 



u 



(57) 



where we neglected terms suppressed by powers of either ^ or p. The Borel function 
which generates all the log-enhanced terms in the large (3q limit is therefore given by 



B{u, x = 1 - e, p) = ^ |- (1 - tx) (p/4)-" r 



■2u 



sm T\u 
2ttu 



1 1 

+ 



1 — M 2 — u 



(58) 



Finally we perform the x integration and incorporate virtual corrections through the 
subtraction of the result at = 1. In spite of the fact that was derived assuming 
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the hierarchy p ^ ^ ^ 1, it is sufficient to get the exact result for the Sudakov exponent, 
since the region of ^ ^ p yields contributions which are suppressed by powers of p. The 
resulting Borel representation of the Sudakov exponent is 



with 



\'D.a{N, g^, 



B^{u,N,m'/q') = e' 



Po Jo 



(59) 



' ■ (1 - u) r(-2u 

. u 



sinvru /2 



2nu 



4.2 Factorization 



Our result (jKHjl for the Sudakov exponent of the e"'"e^ — > Q + X differential cross section 
involves two external scales, q^ and m?. Therefore, in addition to Sudakov logs L = In A^, 
it contains collinear logs / = hiq^ /m?. The latter can be seen a reflection of evolution 
from the scale at which the fragmentation process of the heavy quark takes place 
(Sec. 2 and 3) to the e+e~ centre-of-mass energy squared q^ at which the hard interaction 
takes place. The contributions to the cross section from these two subprocesses can be 
written separately by introducing some factorization procedure. Since Sudakov logs 
emerge from both subprocesses, this separation will be useful to distinguish between 
them. 

Factorization takes the form of a product in moment space, so the Sudakov exponent 
can be written as a sum of the following three exponents. 



ln(T(A^, g^, 



DGE 



In J{N, q^ 



DGE 



+ \nE{N,q^,m^ 



+ lnL)(Ar,m^ 



DGE 



(61) 



where \n.D{N,iv?) corresponds to Eq. (03}, i.e. the (process independent) heavy-quark 
fragmentation depending on the scale m^. In J{N, q^) corresponds to the e^e~ coefficient 
function depending on the scale g^, and lni^(A^, g^, rn?) represents the evolution between 
these two scales and thus it depends on both. 

Performing factorization, the Sudakov exponent corresponding to the fragmentation 
subprocess is the evolution is given by 



\\iE{N, g^m^ 



= 9l 

DGE /3q 7o 



du T{u) (gVA') 



Ba{u) 
u 



IniV 1 



and the Sudakov exponent of the e^e coefficient function is 



lnJ(iV,g2)| = -^^ r duT{u) (<f / 

DGE Po Jo ^ ^ 



(62) 



(63) 



sinvru /2 



- + 



2'Ku \u 1 — u 2 — u 



u 
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The essential ingredient in obtaining a factorized formula is that each of the three ex- 
ponents corresponding to the separate subprocesses has a well-defined perturbative ex- 
pansion: there are no 1/u singularities in these Borel functions. Of course, there is some 
arbitrariness in this procedure as far as higher-order terms are concerned. Here we chose 
the natural factorization scales as the external scales, and q^. The subtracted term 
Ba{'^) = 1 + O'2'u + . . . depends on the factorization scheme. In MS, the coefficient 02 
is given by (021) • Owing to the fact that the expressions above correspond to all-order 
resummation, factorization-scale and scheme dependent terms cancel out completely in 
the sum (jfill) . Of course, if any of the three exponents is replaced by some fixed or- 
der or a fixed-logarithmic accuracy approximation, some factorization-scale and scheme 
dependence will appear in the final result for the cross section. 

Note that while J(iV, q^) and D{N, rn?) contain Sudakov double logs (upon expanding 
the exponent, the highest power of In is twice that of the coupling at each order), the 
evolution factor E{N ^m?) contains at most InA^ to the same power as the coupling. 
This is so because lnii^(A^, g^, m?) has just one power of InN to any order in |Hn|l4()pi-{5j . 

Having performed factorization we can compare the results obtained for the two 
subprocesses fl44|) and ()63|) to the corresponding process-independent calculations. The 
DGE Sudakov exponent for the heavy-quark fragmentation subprocess has been identi- 
fied with the result obtained in the process-independent calculation (jUj). It follows that 
the only scale which plays a role in the fragmentation process D{N,m?) at large is 
m/N . As shown in Sec. 2, this is a general property of the fragmentation function which 
holds beyond the perturbative level. 

Similarly, the DGE Sudakov exponent for the e"'"e~ coefficient function (j63|) can be 
identified as the one appearing in a massless quark fragmentation process which was 
calculated in [22j ■ The same all-order result was obtained p H l23 | l2()j for the F2 structure 
function at large A^. As discussed in P^ l^. In J(A^, g^) ^ is associated with radiation 
from the undetected quark (or the outgoing quark in a deep-inelastic process at large 
Xbj) and the formation of a jet (this explains the notation J) with invariant mass q'^/N, 
which recoils against the detected quark. This identification means, in particular, that 
In J{N, q^) is not specific to the particular process in which two heavy quarks are 
produced: it would be the same if the recoiling jet would have been formed around a 
light quark. The internal structure of the jet is not resolved as the relevant scale is its 
invariant mass. 

The fully factorized process of heavy-quark production in e"'"e~ annihilation is shown 
in Fig. The hard blob represents radiative corrections with virtualities of order of 
the centre-of-mass energy squared q^. The detected heavy meson (pi) emerges from 
the lower jet in the figure. Although a priori two-jet production is symmetric, the fact 
that one measures a single particle inclusive cross section breaks the symmetry: the 
measurement is sensitive to different features of the two jets. The fragmentation blob 
represents D{N,im?). It accounts for soft and coUinear radiation on the scale m/N . The 
recoiling jet on the upper part of the figure represents J{N, g^), accounting for soft and 
collinear radiation of the scale q'^/N. The soft blob is associated with the evolution factor 
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E{N,q'^ ,m?). It accounts for gluons whose momentum components are all small. These 
gluons interact with eikonalized quarks: They are insensitive to the invariant mass or 
the structure of either of the two jets. 



Jet 



T^(q) 




Fragmentation function 



Figure 5: Factorization of the process 7* — ^ Q + X aX large A^. 



Alongside the dominant perturbative corrections which are summarized by ()6H). the 
large- limit singles out certain non-perturbative corrections. Such corrections appear 
in two of the subprocesses described above: in D{N,m'^) they appear on the scale m/N 
and in J{N^q'^) on the scale q^/N. The evolution factor E{N,q'^,m'^), on the other 
hand, has no power corrections. The presence of power corrections, as well as their 
parametric form, can be deduced from the renormalon singularities in the corresponding 
Borel sums. In both D{N, m^) and J{N, g^), the dominant corrections at large N appear 
in the Sudakov exponents, Eq. and (j63|) . so these power corrections exponentiate 
together with the perturbative logs. The renormalon singularities of In J{N, q^) 
were analyzed in [2111231201 • They lead to h?N/q^ and A^A^^/g^ corrections at the 
exponent, as summarized by Eq. (50) in [221. As discussed in the previous section the 
renormalon singularities of the heavy-quark fragmentation exponent \n.D{N,m?) 
appear at u = 1/2 and at any integer and half integer such that u > 3/2. There is no 
renormalon pole at m = 1 owing to the factor {1 — u) in ()45|) . 

5 Implications for phenomenology 

In the previous sections we studied the problem of heavy-quark fragmentation, con- 
centrating on the large- A^ limit. We began by identifying the formal limit in which 
simplification occurs, namely where A^ and m get large simultaneously with the ratio 
m/N fixed. We then performed a process- independent perturbative calculation by DGE, 
facilitating the resummation of Sudakov logs as well as that of running-coupling effects. 
The combined treatment of Sudakov logs and renormalons sets the basis for a systematic 
parametrization of power-suppressed corrections on the scale m/N, which exponentiate 
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together with the perturbative logs. Then, speciahzing to the case of e^e~ annihilation, 
we demonstrated how this result for the Sudakov exponent of the fragmentation func- 
tion emerges from the process-specific renormalon calculation upon factorization. The 
Sudakov-resummed coefficient function in the e^e~ case has been identified with the 
familiar jet function which plays an important role in light-quark fragmentation and in 
deep inelastic structure functions at large x. 

The purpose of this section is to study the implications of these new results on phe- 
nomenology. While there are many possible applications, we concentrate here on the one 
where most precise data is available, namely bottom production is e^e~ annihilation at 
LEPl. In this study we work directly with data in moment space. This is a natural 
choice from a theoretical perspective, but, as we will see, not an easy one from the ex- 
perimental point of view: the presence of strong correlations between different moments 
requires a very careful treatment of the errors. 

In our phenomenological analysis we will concentrate on the large- a: region. This 
region is inaccessible with a perturbative approach without Sudakov resummation. It 
is still very problematic when standard NLL resummation is applied JHl; as the per- 
turbative spectrum is significantly more peaked than the data and it becomes negative 
near x = 1. In such circumstances, matching the perturbative result with the non- 
perturbative contribution becomes awkward. 

This is not intended to be an exhaustive phenomenological analysis of the heavy- 
quark fragmentation function. We concentrate on improving the resummation at large x 
and the parametrization of power corrections on the scale m/N, and we simplify other 
aspects. We do not take into account power-suppressed corrections in m'^/q'^, nor do we 
include O(a^) terms which are free of Inm^/g^ and IniV enhancement. Such effects have 
been considered elsewhere j2Z], and found to have a limited impact. 

In Sec. 15.11 we study the implication of our approach at the perturbative level, com- 
paring it to NLL Sudakov resummation. We also detail there the prescription we use for 
matching the resummed result with the fixed-order calculation and the way we deal with 
infrared renormalon ambiguities. In Sec. 15.21 we address the non-perturbative contribu- 
tion to the fragmentation function. We first discuss the data and the way we use them, 
and then present results of various fits where and the non-perturbative parameters 
are determined by the data. 

5.1 The perturbative result 

5.1.1 Factorization, matching and regularization of renormalon singularities 

Making full use of the concept of factorization as outlined in |T3] , we write the normalized 
Mellin moments of the e~^e~ differential cross section {1/ a){d(T / (Ixe) as a product 

a(iV, q\ m") = C{N, q'; E{N, fx% /i^^) D(iV, m'; /i^^) , (64) 
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where D{N,m'^; fi^p) is the process-independent heavy-quark fragmentation function'' 
defined in Eqs. ([T)) and (jH), C{N,q'^, fij^) is an e~^e~ massless coefficient function, and 
E{N, fip, flip) is an MS Aharelh-Parisi evolution factor, given in Eq. (43) of [THj, which 
is associated with the ultraviolet singularity of D{N,m?] fi^p) (see e.g. and the 

collinear singularity of C{N,q^]fLy). We will choose the natural factorization scales: 



^4 



and (i^p 



m 



Varying these scales was shown to have a small effect. 



provided that resummation is performed in each of these functions. 

In contrast with the evolution factor, the fragmentation function D{N,m'^) and the 
coefficient function C{N,q'^) contain double logs of N and infrared renormalons. If 
q'^/N <^ the coefficient function can be safely calculated using Sudakov resummation 
to NLL accuracy, as done in jT^]. The same holds for the fragmentation function if 
m/N <C A. The latter, however, holds for the first few moments at best. Therefore, 
here the effects of renormalon resummation and the corresponding power corrections are 
expected to be important for phenomenological applications. 

The Sudakov exponents corresponding to the fragmentation function and the co- 
efficient function are given to all orders in the large-/3o limit in Eqs. ()44|) and ()63|). 
respectively. These exponents can be employed in an expanded form, e.g. truncated at 
NLL order. 



In D{N, 



q^^^- In + q-^^.^ 

ymi ' I i/mi 



(65) 



where g>^^ and g^^-^ are given in eqs. (74) and (75) of Ref. jTHj, or evaluated in full, thus 
providing some resummation of subleading logs. In either case, we match^ the exponents 
to the fixed-order results at 0{a^ by the so-called 'log-i? matching' procedure: 



DGE 



fixed order 



lnD(Ar, m^) 



overlap 



(66) 



where In i) ( A^, m} 



DGE 



is given by 



lnD(Ar, 



fixed order 



is the expansion of the 



logarithm of the full result^ for D{N,'m'^) truncated at 0{as) and 



\nD{N, 



overlap 



- In^ A^ + (1 - 27s) In AT 



(67) 



is subtracted to account for the terms which are double counted when adding the two 
previous contributions. Note that we suppress here the explicit dependence on the renor- 
malization and the factorization scales (see Eq. (45) in ^^), which are both set equal 
to m. The matching of the coefficient function C{N, q^) is done in the same way: 



\nC{N,q^) = In J {N,q^ 



+ In C{N,q'^ 



fixed order 



- In J{N, q^ 



overlap ' 



(68) 



'^D{N, m^; fiQp) is sometimes called the "initial condition" for the evolution. This term was therefore 
given the label "ini" in |13j . 

^Note that the matching procedure used in I :'> is somewhat different. The numerical differences are, 
however, negligible. 

^This result was first obtained in |H]; see Eq. (A. 13) there. It can also be read off Eq. (|43|l of the 
present paper. 
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where In J(A^, g^) ^ is given in (jUHj) and C{N,q^] n^), which includes contributions to 
the hard subprocess as well as the jet (see Fig. E)) is given at 0{as) in Eq. (A. 12) of |Hj. 

As explained above the standard factorization procedure is associated with poles at 
M = in the Borel representation of the separate subprocesses. For the fragmentation 
function this is a logarithmic ultraviolet divergence and for the coefficient function a log- 
arithmic divergence of collinear origin. With the appropriate subtraction both functions 
become free of u — > singularities. Thus, each of them has a well-defined perturba- 
tive expansion (which is factorization scheme dependent) to any order. However, this 
expansion does not converge owing to the factorial increase of the coefficients induced 
by renormalons: In the Borel representation infrared renormalons show up as poles at 
positive u, which need to be avoided when performing the integral. This leads to a 
power-suppressed ambiguity. 

Let us examine the nature of the perturbative expansion of the Sudakov exponents. 
The DGE Sudakov exponent for the fragmentation function \nD{N ,m^) of Eq. (03)), 
expanded in powers of A{m?), is evaluated in Fig. IHl (lower right box) for m = 4.75 GeV, 
Q,Ms (jYi'i'^ = 0.218 and = 4, 12 and 30. The plot shows progressive partial sums of a 
fixed-logarithmic accuracy, starting at LL, (m=0 in the figure) and proceeding to higher 
logarithmic accuracy. It should be noted that the result is exact only to NLL, while at 
higher orders the values correspond to the DGE extrapolation from the large-/3o limit. 
The renormalon effect leading to divergence of the series sets in early on: the minimal 
term is NLL (m=l) in all three cases. 

Similar plots of the Sudakov exponents In J{N, g^). In E{N, g^, m^) and lna(A^, g^, m^), 
of Eqs. (jfiHj) ■ (jn21) and (fHn|) . respectively, are shown in the other boxes for a center-of-mass 
energy of \/q^ = Mz and af^ (^I) = 0.118. Here the fixed-logarithmic accuracy expan- 
sion is in powers of A{M'^). It is clear that the divergence of the coefficient function is 
rather mild: In J does not reach the minimal terms up to m=8; hiE converges and the 
Sudakov exponent of the total cross section In a reaches the minimal term at NNLL or at 
N^LL depending on A^. Here the divergence is induced by the fragmentation subprocess. 

When the expansion is truncated, e.g. at NLL order, the renormalon ambiguity does 
not appear but power accuracy is usually not reached. If one accepts that the DGE 
result represents well the contribution of subleading logs, going to power accuracy 
simply requires avoiding any arbitrary fixed-logarithmic accuracy truncation. Since the 
series does not converge, a sensible possibility is to define the perturbative sum by 
truncation at the minimal term, i.e. just before the series starts diverging. The minimal 
term scales as a power. This definition was recently used ^.20^ in the context of DGE for 
the deep inelastic structure function F2, where a purely perturbative result was shown 
to be consistent with the data. However, in the case of the heavy-quark fragmentation 
considered here power corrections are particularly large and the series starts diverging 
already at low orders. In this case truncation at the minimal term would be far from 
optimal as it leads to discontinuities in the perturbative result as a function of A^ whose 
magnitude is comparable to the power correction itself, thus inducing 0{1) discontinuities 
in the non-perturbative contribution. The divergence of the coefficient function C{N, g^) 
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Figure 6: Evaluation of the Sudakov exponents with increasing logarithmic accuracy. 
The minimal term in the expansion is denoted by a cross. 



is much softer, so here truncation at the minimal term would be numerically sensible. 

A better regularization prescription is to evaluate the Borel integral in Eq. ()44j] 
directly, with a Principal Value (PV) prescription applied at each pole. Fortran routines 
exist for evaluating such integrals efficiently provided that the exact position of the 
singularity is known, as in our case. This will therefore be our default prescription. All 
the numerical results for DGE we present below are obtained in this way. 

5.1.2 Comparison between DGE and NLL resummation 

We now compare the numerical results for NLL and DGE resummation with PV regu- 
larization. Our default parameters will be = 91.2 GeV, m = 4.75 GeV, Nf = 5 and 
A = 0.243 GeV, corresponding to af^ (M|) = 0.118 and af^ (m^) = 0.218. We will be 
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using a two-loop running coupling^° given by eqs. (jH^ and (jHHj) . or, explicitly, by 



/So V 1 + 



where 



w{s) = W-i (- exp(-s/5 - 1)) , 
Here W stands for the Lambert W function |l2l 
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Figure 7: Left: the matched heavy-quark fragmentation function D(A^, m^) with DGE 
(solid) and NLL (dashes) resummation. Right: the same for the e^e~ coefficient function 



Figure [3 shows a comparison between the NLL and DGE resummed results for the 
heavy-quark fragmentation function D{N ^vp?\ and for the e^e~ coefficient function 
Ci^N^q^) (both are matched according to and to the full 0{as) result). One 
can immediately see that the NLL result for D^N^m?) breaks down at = = 
exp(7r/ (2/3oas(m^)) ~ 40 where the 5^;^] functions become singular. The DGE result with 
PV regularization, on the other hand, remains moderate and shows no singular behavior. 
Indeed, owing to the procedure adopted where the Borel integrals are directly evaluated 
rather than expanded the DGE result is free^^ of Landau singularities. 

This does not mean, of course, that the perturbative DGE result by itself is physically 
meaningful: power corrections are needed. The advantage is, however, that this result, 
being free of spurious singularities, provides a good basis for the parametrization of 

^°Note that the exact two-loop evolution Eq. is close but not identical to the commonly used 
expansion of the running coupling in In/i^, (see e.g. Eq. (29) in The choice we have made for 

obtaining our default value a^(M|) = 0.118 corresponds to A = 0.226 GeV when Eq. (29) in ^ is 
used. 

^^Note that the convergence of H44|l at u — > 00 for any N is guaranteed by the factor r(— u). This 
factor also introduces an infinite set of renormalon singularities. 
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power corrections for any A^. In the case of C{N, q^) the difference between the NLL and 
DGE results is much smaller than for D{N, m^). The Landau singularity in the Sudakov- 
resummed coefficient function is located much further, at = exp(7r/(/5oQ;s(g^)) ~ 10^. 
Thus here subleading logs are smaller and power corrections can definitely be neglected. 




Figure 8: Left: a comparison between DGE and NLL for heavy-quark production in the 
e^e~ process. Right: the same curves m. xe space. ALEPH jl^ data are also shown. 



Figure IHl shows a similar comparison for the full e^e~ cross section a{N,q'^,m'^) of 
Eq. (jMj) . For reference we show also a curve with no soft gluon resummation as well as 
the experimental data from ALEPH The theoretical curves in space are obtained 
via numerical Mellin inversion: 

D{x) = — f dN x-^b 
2Tii Jc 

In the NLL case the Minimal Prescription ^1] is used: the contour C (see Fig. ^ passes 
between the origin and the leftmost Landau pole. In the case of the (PV-regularized) 
DGE result a specific prescription is not required since there are no Landau singularities. 
Any contour passing to the right of the origin would yield the same result. This may be 
surprising at first sight, since the original expression for the fragmentation function in 
X space (j30|l does have a convergence constraint (see for example p6|l ). However, as we 
saw, when going to moment space this constraint is traded for renormalon singularities. 
Landau singularities instead appear only if the resummed expression is expanded in 
terms of the coupling. The inherent ambiguity of the perturbative description of the 
fragmentation function as x — > 1 appears in the Borel sum as a renormalon ambiguity: 
In order to compensate for it one needs to specify an infinite set of power terms (A^A/m)". 

A hint about the size of such contributions can be obtained by considering different 
prescriptions for performing the Borel integration, PV being one of the possible choices. 
Evaluating lnD(A^, m^) in (jl^ by going above and below the first pole at m = 1/2 we 
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Figure 9: A prototype of the contour C employed in the numerical Mellin inversion to 
X space. 



obtain the band bounded by the dot-dashed curves in Fig. |H1 The experimental data turn 
out to be within the band: This supports our expectation that one can describe them 
by modifying the resummed perturbative result with the first few power corrections. 

5.2 Non-perturbative contributions 

In general, the non-perturbative contribution to the fragmentation function should be 
included on top of the perturbative one through a convolution in x. Moments of the B- 
meson cross section as measured at LEP can then be written as the product of the pertur- 
batively calculated moments a^'^ {N , , m?) and the moments of some non-perturbative 
function. 

Our first result is that at large and m this function depends on the ratio NA/m 
alone, up to corrections of order Unfortunately, with just one heavy flavour (charm 

is probably too light) the predictive power of this statement is largely lost: fixing m the 
fragmentation can anyway depend only on A^. This situation can be contrasted, for 
example, with the case of event-shape distributions where non-perturbative corrections 
can a priori depend on the centre-of-mass energy as well as on the shape variable: Thus 
the statement that the corrections can be described by a shape function of a single 
argument is already quite constraining PH ll9j. 

The emphasis in the phenomenological analysis is therefore on the particular depen- 
dence of the fragmentation function on NA/m which is predicted by the renormalon 
model (jlHjl . Let us write 



a{N, q\ M') ^ a^-{N, q\ m") DfMN - l)A/m) 



(72) 



where 




(A^- 1)A 



) 



n 



m 



(73) 
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The separation being, of course, ambiguous, the perturbative part in Eq. (f?^ describes 
the production of a bottom quark surrounded by a cloud of soft gluons and the non- 
perturbative one describes its hadronization into an observable B meson. As before, we 
have not explicitly shown the dependence on the QCD scale A which the perturbative 
part o"""^ acquires via the strong coupling. The dependence of a{N, g^, M^) on the heavy 
meson mass M rather than the heavy quark mass m serves as a reminder that it refers 
to the observed particle. 

In Eq. (|73|) we have made a few modifications compared to Eq. ()48|1 . Firstly, the 
factor Cp/Po has been included into the parameters e„. Secondly, we have introduced a 
term proportional to (A^A/m)^, which is absent in the renormalon result. By including 
it and fitting to data we can examine this feature of the renormalon model ()48|) . Finally, 
we have replaced A^ — *• A^ — 1. This modification is of course allowed in the large- 
A^ limit, and it implements the constraint that for A^ = 1 (i.e. the normalized total 
cross section) the non-perturbative correction vanishes. Eq. (f7^ will be the master 
equation for our phenomenological analyses of the LEP data. The default choice for the 
perturbative calculation of the 6-quark differential cross section a^'^{N ,m?) will be 
the PV-regularized DGE. 

5.2.1 The experimental data 

Several experimental collaborations have recently published high-statistics and high- 
accuracy data for i?-meson production in e^e~ collisions at the Z° peak, i.e. centre- 
of-mass energy around 91 GeV. Most of the data have been published in the form of 
differential distributions in the variable xe-, representing the ratio between the i?-meson 
energy in the laboratory frame and the beam energy. The formal definition of the frag- 
mentation function (0) is in terms of the longitudinal momentum fraction x. However, 
for large the two coincide: — up to corrections (jTfjl of order vr? which we 
neglect. From here on the two notations will be used interchangeably. 

Data in x space, along with their covariance matrices, have been published by the 
SLD gni, ALEPH Hg, DELPHI \^ and OPAL [TT] Collaborations. The DELPHI 
Collaboration has also published moment data, including the very important correlations 
between the moments up to A^ = 6. 

Rather than converting our moment-space expressions to x space as usually done, we 
shall perform our fits directly in moment space, using the CERN Library minimization 
routine MINUIT |1H]- Given the x-space data and covariance matrices as inputs, moments 
and their own covariance matrices can be calculated. The integral (jlj) defining Mellin 
moments can be replaced by a discrete sum over all the bins in x space (ntina in total), 

a(A^,g^M2)^ ^xf-^Ax./,, (74) 

where /j is the normalized value of the i-th bin, Xi its central abscissa and Axj its width. 
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This equation represents a linear transformation of the vector /j, 

£ = A/, (75) 

and the matrix A is defined by Aj^ji = xf'^Axi. Given the covariance matrix F for the 
fi values, one can build the covariance matrix S for the cr{N, q^, M^) moments as 

S = A F . (76) 



N 


ALEPH 


DELPHI 46,49 


DELPHI 


2 


0.7163 ± 0.0085 


0.71422 ± 0.0052 


0.7147 ± 0.0045 


3 


0.5433 ± 0.0097 


0.5401 ± 0.0064 


0.5413 ± 0.0057 


4 


0.4269 ± 0.0098 


0.4236 ± 0.0065 


0.4248 ± 0.0060 


5 


0.3437 ± 0.0096 


0.3406 ± 0.0064 


0.3419 ± 0.0059 


6 


0.2819 ± 0.0094 


0.2789 ± 0.0061 


0.2804 ± 0.0057 


7 


0.2345 




0.0091 




0.2333 




0.0054 


8 


0.1975 




0.0087 




0.1965 




0.0050 


9 


0.1680 




0.0084 




0.1672 




0.0047 


10 


0.1441 




0.0080 




0.1435 




0.0044 


11 


0.1245 




0.0076 




0.1241 




0.0041 


12 


0.1084 




0.0072 




0.1081 




0.0038 


13 


0.0949 




0.0069 




0.0947 




0.0036 


14 


0.0835 




0.0065 




0.0835 




0.0033 


15 


0.0738 




0.0062 




0.0740 




0.0031 


16 


0.0656 




0.0058 




0.0659 




0.0029 


17 


0.0585 




0.0055 




0.0589 




0.0027 


18 


0.0524 




0.0052 




0.0529 




0.0026 


19 


0.0471 




0.0049 




0.0477 




0.0024 


20 


0.0425 




0.0047 




0.0432 




0.0023 


21 


0.0384 




0.0044 




0.0393 




0.0022 



Table 1: Experimental data for normalized moments of weakly-decaying S-meson pro- 
duction in e~^e~ collisions from the ALEPH and DELPHI Collaborations at LEP, eval- 
uated according to Eq. (f7^ (left and right columns). The set in the central column 
corresponds to the five moments published by DELPHI in jlH] (the N = 2 point does 
not include here the correction for initial state electromagnetic radiation |12j). 

Table C] contains most of the moment-space data which will be used in the fits. 
Moments constructed from ALEPH and DELPHI data are shown. Moments higher than 
those shown in Tabled can of course be calculated, and they have also been used in some 
of the fits. The amount of independent information they carry is however limited. 

The central values reported in Table Q show a very good agreement between the 
ALEPH and DELPHI sets. DELPHI moments are reported twice in the Table. The 
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data in the central column are directly provided by jlHllini, where they were pubhshed 
in preliminary form, whereas we have calculated the ones in the right one according to 
Eq. ()74|1 . It should be noted that the DELPHI Collaboration extracted the moments 
directly from the unfolding program they used, rather than calculating them by summing 
up the published bins in x space. This accounts for the small differences. The two sets 
of moments are, however, compatible within errors. For full consistency, we shall make 
use of the published set of moments with the published covariance matrix, and of our 
own values together with a covariance matrix we extract ourselves according to Eq. (|76|l . 
The latter also presents some differences with the published one. 

As a test of the fitting program, we have reproduced the fit that the DELPHI Collab- 
oration performed in jlH] to the five correlated moments which were published there, em- 
ploying a so-called Kartvehshvili [H] functional form, D^^{a^x) = {a + l){a + 2)x°'{l—x). 
Using their parameters, their moments and covariance matrix, and NLL resummation 
only, we do recover their result: a = 17.07±0.57, = 114.7 for 5-1 de grees of freedom. 
However, fitting the moments we calculate ourselves (right column in Table and using 
our own covariance matrix, we get instead a = 14.30 ± 0.37 and = 159.8. The appar- 
ent incompatibility of these results should mainly be attributed to the intrinsic failure 
of the fit, as indicated by the very large values. It also shows, however, how strong 
is the effect of small variations in the moments and in their covariances, due to the high 
degree of correlation and to the very small experimental errors. We have also performed 
the same fit to the first five ALEPH moments, obtaining a = 21.5 ± 2.0 and = 15.1. 

We wish to point out that the very strong correlations between the experimental 
data for the moments are an important feature, which cannot be neglected when fitting 
a hadronization model. Extracting a meaningful covariance matrix for many moments 
at once is however no straightforward task. Because of the strong correlation between 
successive moments, inversion of such a matrix - needed when evaluating the fo^' 
the fits - quickly becomes an intractable problem. In practice the numerical accuracy 
issue becomes acute when more than five or six consecutive moments are fitted together. 
Possible solutions are to use only a few moments or to pick non-consecutive moments. 
By means of these choices the matrix becomes smaller and the degree of correlation is 
lessened. One observes that matrices containing five or six moments, taken two or three 
moments apart, can be inverted with acceptable precision up to ~ 20 . Using only 
three or four moments a sufficient numerical accuracy can be obtained up to A^ ~ 30. 
We make use of this expedient in some of the fits described below. 

5.2.2 Fits to low moments 

In order to test the renormalon model for the non-perturbative fragmentation function 
we would like to fit at least three parameters: A, ei and €2- ei is the leading power 
correction corresponding to the shift of the entire distribution in x space in units of 
A/m. It is expected to be of order 1. Based on the renormalon model ()48j) . £2 is 
expected to vanish. Note that varying A, i.e. the value of af^ (M^), influences both the 
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perturbative and the non-perturbative parts of the cross section ((721). approach to 
heavy-quark fragmentation can thus be tested by checking that the fit returns a value 
which is consistent with other determinations of af^ (^^1)- 

With this task in mind and the data described above at hand we should choose a 
subset of moments to be fitted. Theoretical limitations exist both for very small and 
for very large N. In the former case neither the soft-gluon perturbative resummation nor 
the power corrections of (jlH|) are expected to dominate. In the latter, parametrizing just 
a few power corrections in the exponent may not be sufficient, as NA/m is not small. 

In this section we concentrate on moments up to = 6. In the next section we 
extend the analysis to higher moments. We use the correlated moments obtained from 
the ALEPH data and shown in Tabled DELPHI moments have also been used for cross 
checks and comparisons but, due to their preliminary nature, we refrain from quoting 
these results. 





e„ = for n > 3 




1.23 ± 0.23 


1.03 ± 0.13 


2.9 ± 1.5 




(fixed) 


0.011 ± 0.095 


-1.1 ± 1.2 


A (GeV) 


0.217 ± 0.025 


0.243 (fixed) 


0.138 ± 0.045 


xVd.o.f 


6.8/3 


7.9/3 


3.7/2 



Table 2: Results of fits to ALEPH moments from = 2 to = 6. 

First, we wish to extract A and test the prediction that the leading power correction 
is of the order of A/m. We do this by fitting A and ei, setting = for any n > 2 in 
Eq. ()73|1 . The results of this fit are shown in the first column of Table |2l The of the fit 
is reasonable albeit somewhat large. It is important to notice that the central value for 
A, corresponding to af^ {M^) — 0.117, is fully compatible with present determinations 
of the strong coupling. The best fit value for ei is of order 1, as expected. The resulting 
curves and the ALEPH data are shown in Fig. ^[ both in A^ and xe space. The overall 
description of the data is reasonable. This is quite remarkable given that, in xe space, 
the theoretical curve is no more than a shift of the PV-regularized DGE perturbative 
result. 

The self-consistency of this fit can be probed by fitting now €2 while fixing A = 
0.243 GeV. Being this value for A compatible with the one returned by the previous fit, 
we expect to get a compatible value for ei. The results are shown in the second column 
of Table El ei is indeed consistent with the one shown in the first column, the is 
similar, and 62 turns out to be consistent with zero. 

Finally, we attempt fitting three parameters, ei, €2 and A, simultaneously. The results 
are shown in the third column of Table El Consistency with previous fits appears rough, 
and the errors are extremely large. 

Because of the problems encountered in the three-parameter fit, and since we ex- 
pect the perturbative corrections we resum and the power correction we parametrize to 
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Figure 10: Left: results of a fit for A and ei (setting e„ = for n > 2) to ALEPH 
moments = 2 to 6. Right: the corresponding curve in xe space, compared to the 
ALEPH and DELPHI data. 

become dominant at larger A^, we now proceed to analyze higher moments. 
5.2.3 Fits to high moments 

In this section we wish to explore the possibility of fitting Eq. ()72|) to large- A^ correlated 
moments. The large-A^ region is challenging in several respects. Considering the data, 
the finite binning in X£;-space limits the amount of independent information contained 
in large- A^ moments. On the theoretical side, as soon as the condition NA/m -C 1 is 
violated the non-perturbative contribution is expected to become comparable to, and 
eventually to dominate, the perturbative one. Moreover, the exponentiation of power 
corrections can no longer be ignored and the specific form of (jUI), e.g. the vanishing 
of the second power of NA/m, becomes relevant. This makes this exercise particularly 
interesting: The validity of our model can really be tested. Yet, it is a priori not known 
up to what value of A^ Eq. (jIH|) with just a few parameters might hold: One should 
keep in mind that when going to extremely large A^ the number of relevant parameters, 
corresponding to increasing powers of NA/m in the exponent, increases, and the fitting 
procedure may get out of control. 

The moments we fit here are those we reconstructed from the X£;-space data and 
covariance matrices^^ published by ALEPH. The fine binning in xe space of the ALEPH 
data helps in providing a good description of the large-A^ moments in which we are 

^^We use covariance matrices with more significant figures [101 than pubhshed in j48| . 
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6,8,10,12,14,16 


13,15,17,19 


7,10,13,16,19 


21,24,27, 30 


ei 

A 

1 A. 

xVd.o.f 


0.831 ± 0.041 
-0.037 ± 0.007 

261 -1- 01^ 
11/3 


0.77 ± 0.15 
-0.024 ±0.028 
955 +0 025 
8.3/1 


0.96 ± 0.12 
-0.063 ± 0.025 
n 2S8 + 022 
9.0/2 


0.817 ± 0.056 
-0.040 ± 0.008 
256 + 020 
3.3/1 



Table 3: Fits to correlated ALEPH moments in the large- A^ region. Only a few non- 
consecutive moments are used, as indicated in the various columns. 

interested^'^. Provided that the correlations are properly taken into account, it allows 
one to consider moments as high as A^ ~ 30. 

In order to achieve a sufficiently^^ accurate inversion of the covariance matrix we 
now make use of the expedient of using non-consecutive moments. The results, shown 
in Table El are remarkably consistent: Over a wide range of moments the best-fit values 
for A, ei and 62 appear very similar, and in line with our expectations. 

We further explore the large- A^ region by employing 'rolling window' fits, using again 
non- consecutive moments. FigurelTTlshows ffis to A, ei and €2, performed within windows 
bordered by iV ± 3, where only four correlated moments are used. At small A^ the 
fits are inconclusive. As we saw above (Table ^ this is a consequence of trying to fit 
too many parameters. In the large-A^ region the situation improves dramatically: €2 
tends to zero, as predicted by the renormalon model, and the value for A stabilizes on 
- 0.25 ± 0.04 GeV, corresponding to af^ (M|) ~ 0.118 ± 0.003, which is compatible 
with other determinations. We stress that this should not be regarded as a reliable 
determination of the strong coupling as the analysis of experimental and theoretical 
errors was not performed. Note the fluctuations in the error bars at large N values, due 
to the deteriorating numerical accuracy. 

For comparison. Fig. ^2 also shows the same kind of fits with NLL resummation. 
The range in A^ is limited here by the non-physical behavior of the resummed result for 
A^ greater than 30 or so, as shown in Fig. |H| The results are clearly at variance with the 
DGE fits: A and ei do not stabilize and €2 does not tend to zero when A^ becomes large. 

5.2.4 Comparison with moment-space and XE-space data 

At the end of the day, one wishes of course to extract a parametrization for non- 
perturbative effects such that, when convoluted with the appropriate perturbative con- 

"'^^For a given binning, at extremely large N the moments are eventually determined only by the very 
last bin at large x, and thus higher moments carry no new information even in an ideal situation where 
the bins are not correlated. This is not yet realized for the ALEPH data at TV of a few tens. For 
example, for = 10 there are twelve bins that contribute more than the last one, for N = 20 this 
number drops to seven and for = 30 there are still five such bins. 

^''Note that the fits in the large-A^ region are performed at the edge of the numerical accuracy 
permitted by the inversion of the covariance matrix. 
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Figure 11: Results of fits to ei, €2 and A (in GeV on the vertical axis) performed with 
correlated ALEPH moments within a "rolling window" , centered at N, where 4 moments 
out of 7 are used. In the left plot the DGE perturbative result (with PV regularization) 
is used. In the right one instead the resummation is truncated at NLL accuracy. 



tribution, it allows for a fair description of the experimental data in both moment and 
xe space. 

The experimental papers [131l43p46|ITF] have tested a number of different parametriza- 
tions, and estabhshed which ones seem to offer a good description of the data. Our 
present attempt differs from theirs in several respects: 

• we only fit data in moment space, and subsequently derive the corresponding xe 
distribution; 

• our perturbative description is based on DGE, matched to the NLO result ()66p. 
and to the NLL Altarelli-Parisi evolution (jMj) . The experiments instead usually 
use the perturbative fragmentation provided by Monte-Carlo programs. We recall 
that a non-perturbative function should only be used in connection with the same 
perturbative description it has been fitted with. Conversely, changing the perturba- 
tive description (or its parameters) can lead to fit vastly different non-perturbative 
functions; 

• we use a renormalon- motivated functional form, Eq. ()73|) . for the non-perturbative 
function, rather than phenomenological models like the commonly used Peterson et 
al. [Zj or Kartvelishvili et al. [B]. Besides the theoretical reasons for using Eq. (|7Hjl . 
we have also checked that the Kartvelishvili or Peterson functional forms do not 
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Figure 12: Left: results of a fit for A, ei and 62 and to ALEPH moments N = 30, 31, 32, 
setting e,i = for n > 3. Right: the corresponding curve in xe space, compared to the 
ALEPH and DELPHI data. 



provide a good description of the data when used with our DGE-improved pertur- 
bative result. 

For illustrative purposes we perform a further fit using exclusively very high moments, 
= 30, 31 and 32. The main disadvantage of such a fit is that higher power corrections 
on the exponent, which are not parametrized but simply set to zero, may in fact play 
some role. On the other hand it has several advantages: 

• one can be confident that the perturbative corrections resummed by ()44|) and the 
power correction parametrized by (jlHj) indeed dominate. 

• the stability observed in our previous fits at large A^ guarantees that the result will 
be independent of the particular set of moment chosen. 

• it ensures a good description of the large- A^ limit, and thus of the region near 
X — > 1. 

The results are shown in Fig. El The fitted parameters (see the figure) are compatible 
with the ones shown in Table El and in Fig. ^2 ^2, although not vanishing within 
the errors, is very small. In order to check the sensitivity to higher power corrections 
in the exponent, we have performed a similar fit (using A^ = 30, 31 and 32) fixing 
62 = and having A, ei and es as free parameters. The result is visually similar to 
Fig. El but the parameters are somewhat different: A = 0.266 ± 0.020, corresponding 
to af^ (^I) — 0.120, and the parameters of the shape function are ei = 0.86 ± 0.11 and 
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es = —0.011 ± 0.003. These differences with respect to Fig. [T21 represent the sensitivity 
to higher power corrections. 

Despite having used large- data only in the fit, the low-A^ moments are fairly well de- 
scribed by the resulting curve, typically deviating by no more than a few percent^^. One 
can moreover notice that the whole curve in x^; space is well reproduced, notwithstanding 
the fact that it was never directly fitted: Its shape is produced by the DGE-improved 
perturbative calculation, properly shifted and modeled by the fitted non-perturbative 
term (see also Fig. |H1 for a comparison with the purely perturbative result). A better 
description of the lower moments and the region left of the peak is expected upon in- 
cluding corrections of order 0{{m'^ / q^)as)i non- logarithmic 0{a1) corrections j2ZI, and 
power corrections of order A/m not enhanced by A^. 

Note that the curve m. xe space remains positive up to xe very close to one, and 
beyond the last experimental point. This constitutes a marked improvement with respect 
to fits based on fixed-order or NLL Sudakov-resummed perturbative results. In these 
approaches a convoluted non-perturbative function can be fitted to the data, and the 
resulting curve describes well the experimental points in the region below the peak. 
However, the cross section turns out negative over a broader region between the peak 
and = 1, and cannot therefore describe the last few data points. We recall that 
this improvement has been obtained here by reGning the perturbative prediction, adding 
leading and subleading Sudakov logarithms within the DGE formalism and introducing 
an appropriate prescription^^ (PV) to regularize the renormalons. This stands in contrast 
with the approach of [ST] where large logarithmic corrections in the function D{x,im?) 
are simply discarded. 

6 Conclusions 

In this paper we considered the problem of heavy-quark fragmentation in the large-x 
region. We first derived a general relation © between the fragmentation function at 
large A^ and a non-local hadronic matrix element at large light-cone separations. When 
combined with the result on the asymptotic form of this matrix element at large m pij , 
it implies that the simultaneous limit where A^ and m get large together, with the ratio 
between them fixed, has a special status: the asymptotic heavy-quark fragmentation 
function depends on this ratio alone. 

Equipped with this exact asymptotic result we proceeded to evaluate the fragmen- 
tation function perturbatively, resumming large logarithms of A^ as well as running- 

^^It is worth noting that the level of accuracy of a few percent on the low-A^ moments is by far sufficient 
for phenomenologically relevant applications such as the description of heavy-quark hadroproduction as 
performed in j^. 

Although the prescription itself is arbitrary, the ambiguity is of a specific functional form, namely 
that of the power corrections we introduce. Thus, upon performing a fit where these corrections are 
parametrized the ambiguity is fully canceled jl7|: other regularizations simply correspond to a redefini- 
tion of the parameters which control power terms. 
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coupling effect. We first calculated the splitting function of an off-shell gluon of a mas- 
sive quark, from which we derived an all-order result for the fragmentation function in 
the large-/3o limit. This result has some interesting features: in x space it does not con- 
tain any renormalons, but it does have convergence constraints for x — > 1. In moment 
space these constraints arc replaced by an infinite set of infrared renormalons, indicating 
certain non-perturbative corrections to the fragmentation function. The physical inter- 
pretation of this result is clear: the perturbative calculation only gives the probability 
of producing an on-shell quark with an energy fraction x surrounded by a cloud of soft 
gluons, whereas the physical quantity (defined in the full theory) is the probability to 
produce a bound state (a heavy meson) with that energy fraction. The gap between the 
two is filled by these power corrections. 

As expected, at large N renormalon ambiguities depend only on the ratio NA/m. 
Thus, the non-perturbative contribution to the fragmentation function at large x can be 
incorporated through a convolution with shape function of m(l — x). Further details on 
this function can be deduced from the structure of the perturbative result: 

• similarly to the logarithmically enhanced contributions, these power corrections 
exponentiate in moments space; 

• the leading power correction is 0{NK/m). Upon exponentiation it generates a 
shift of the perturbative spectrum in x by an amount proportional to A/m; 

• the second corrections in the exponent, 0{N'^IS? /w?'), is absent, suggesting that 
the shift may be a good approximation in a relatively wide range in x] 

• higher order corrections such as 0{N^h? /m^) which also appear in the exponent 
modify the shape of the spectrum near x ~ 1. 

Regarding the non-perturbative fragmentation function as a set of non-perturbative 
power corrections is useful only if the perturbative sum itself is performed taking Su- 
dakov logs as well as infrared renormalons into account. This is the purpose of DGE. 
The significance of Sudakov logs at large x has long been understood, while that of 
renormalons is more subtle. In applications of Sudakov resummation in the asymptotic 
regime (e.g. in the coefficient function for a single-particle inclusive cross section at large 
- see Sec. 4 and 5) power corrections arc negligible, and a fixed-logarithmic accuracy, 
e.g. NLL, is sufficient. For the heavy-quark fragmentation function, on the other hand, 
power corrections are essential, and therefore so is the proper separation between the 
perturbative sum and the power corrections. A proper separation can be achieved only 
if renormalons, which manifest themselves as subleading Sudakov logs, are resummed. 
Furthermore, we found that dealing with the renormalon ambiguity by means of a PV 
prescription, one avoids Landau singularities and obtains a sensible asymptotic behavior 
at large N . 

The numerical analysis demonstrates the significance of this property. In contrast 
with the fixed-order result or even the NLL resummation, the DGE perturbative result 
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with PV regularization (see Fig. |HJ, does not become large nor negative, but remains 
instead positive up to x = 1. This is a resuh of refining the perturbative calculation for 
heavy-quark fragmentation, by including the effects of leading as well as subleading soft 
logarithms via DGE, while introducing a systematic prescription (PV) to separate the 
perturbative and non-perturbative contributions. Upon convolution with the properly 
fitted non-perturbative contribution, this feature allows for a good description of the 
data up to X ~ 1. 

While concentrating on the large-x region, many other features of the fragmentation 
process have been simplified. We neglected power corrections in m? /q^ as well as pertur- 
bative corrections of order a1 which are not logarithmically enhanced. We also neglected 
power corrections of the type A/m (not enhanced by A^), which certainly play an im- 
portant role at small N . All these aspects definitely deserve further investigation in the 
future. In spite of these simplifications and of the fact the non-perturbative correction 
we employ has minimal flexibility, we obtained a fairly good description (Fig. ITUland lT^ 
of the ALEPH data on 5-meson production. In doing so we extracted a value for as 
which is compatible with other determinations (although a systematic error analysis has 
not been performed here). We also found that the non-perturbative parameters for the 
fragmentation function match the expectations: the shift of the perturbative spectrum 
in X is indeed of order A/m and €2 in Eq. (|7Hjl turns out to be very small. 
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